† Corresponding author. E-mail:
Project supported by the National Natural Science Foundation of China (Grant Nos. 11274168, 11374157, 11174138, 11174139, 11222442, and 81127901) and the National Basic Research Program of China (Grant Nos. 2010CB327803 and 2012CB921504).
In this paper, acoustic scattering from the system comprised of a cloaked object and the multilayer cloak with only one single pair of isotropic media is analyzed with a recursive numerical method. The designed acoustic parameters of the isotropic cloak media are assumed to be single-negative, and the resulting cloak can reduce acoustic scattering from an acoustic sensor while allowing it to receive external information. Several factors that may influence the performance of the cloak, including the number of layers and the acoustic dissipation of the medium are fully analyzed. Furthermore, the possibility of achieving acoustic invisibility with positive acoustic parameters is proposed by searching the optimum value in the parameter space and minimizing the scattering cross-section.
In 2006, Pendry et al. proposed the concept of transformation optics, and used this coordinate transformation technique to design an invisibility cloak for electromagnetic waves.[1] Ever since them, the coordinate transformation technique has attracted worldwide attention. Such a method was extended to the field of acoustics by Cummer and Schurig in 2006.[2] Since then there has been growing interest in acoustic invisibility[3–9] due to the possibility of hiding objects from detecting the acoustic signals, which may lead to applications in a variety of practical situations such as noninvasive acoustic sensors.
Recent research suggested that it is possible to “cancel” the cloaked object as well as to share information from its surroundings by using the complementary objects made of a double-negative medium (DNM).[10] While DNM has been proposed and demonstrated by using acoustic metamaterials,[11–15] the corresponding unit cells are usually related to a resonate effect and consequently are only effective in narrow band. Besides, the resonance-induced dissipation will inevitably harm the performance of the resulting device. Therefore the invisibility cloak based on DNM would be less practical for acoustic waves. Zhu et al. proposed the superlens cloaking by using single-negative medium (SNM) in 2011.[16] Compared with the DNM device, the invisibility cloak based on SNM would generally be easy to fabricate and effective to work in broader bandwidth.[17] Although the superlens cloaking by Zhu et al. is proved to be effective in reducing acoustic scattering from an acoustic sensor, such a device is composed of inhomogeneous anisotropic media which require, even in a reducible way, a large number of isotropic media with different parameters, and thus bringing the difficulty in fabricating it. Later, Zhu et al.[18] revised their folding transformation scheme and simplified the cloaking method, in which only three types of SNMs are needed. Xu et al.[19] further simplified the superlens fabrication process, and in their scheme, only one single pair of isotropic media is required for reducing the acoustic scattering.
In this paper, a recursive numerical method of calculating the acoustic scattering from the structure with a homogeneous isotropic multilayer shell is given in detail, and based on this method the SNM superlens cloaking is fully analyzed. The influences of several factors that may hinder the performance of the cloak, including the number of layers, acoustic dissipation of the media, etc., are investigated. According to the results above, we further investigate the possibility of achieving a positive parameter invisibility device, by searching, in the parameter space, for the optimum value that can render the scattering cross-section minimum.
The rest of the present article is organized as follows. In Section 2, we first introduce the theory of superlens cloaking and then give the theoretical formulation for scattering analyses with the recursive numerical method. In Section 3, we present in detail the numerical results for the cloaking performance in various situations, analyze the factors that may influence the invisibility effect, and give an example of an invisibility cloak based on scattering cancellation. Finally, we draw some conclusions from the present studies in Section 4.
For simplicity, we restrict the problem to a two-dimensional (2D) case. A sketch of the folding space transformation is shown in Figs.
Now we focus on the annular shell B′ which is divided into N bilayers
In order to present the recursive procedure clearly, we begin with the case of only one single cylindrical scattering object. The analysis focuses on acoustic waves propagating along the cross section of a homogeneous infinite cylinder in a host medium. Cylindrical coordinates (r,θ,z) are used throughout the analysis, and the material is assumed to be isotropic. The wave vector is in the plane (r − θ), and the wave-front extends infinitely along the z direction. We assume that the external incident acoustic wave can be described as
We show part of the multilayer structure in the first quadrant in Fig.
The acoustic field can be regarded as a superposition of the standing wave (corresponding to the Bessel function) and the scattered wave outward (corresponding to the Hankel function of the first kind). For instance, in the host medium, it is expressed as
If we already know the incident wave, that is, the value of α0n is given, then there are two general approaches to determining the acoustic scattering field. One is to list all the continuity conditions at all the interfaces, and solve the undetermined coefficient at one time;[5] the other is the recursive approach,[21] which is employed in this paper. The general idea of this recursive method is similar to the transfer matrix method.[22,23] Consider that the acoustic wave in layer j − 1 impinges onto layer j (where j = 1,..., N), the incident, reflected and transmitted waves can be written as
On the other hand, when the wave travels from layer j towards layer j − 1 (where j = 1,...,N), the incident, reflected and transmitted waves can be written as
When only the innermost layer exists, there are two media in the space, where the scattering object is medium N with radius aN and the host material is medium N − 1. When the incident wave is given as
The total reflection coefficient in layer N is represented by
Supposing that the incident wave has the form of
According to the method above, the recursion moves outward in the descending direction of j until the outermost host medium becomes an actual host medium, and then we can obtain the whole distribution of acoustic wave field of the multilayered scattering structure. The details of this solution procedure are as follows.
First, from the innermost layer, the multiple scattering factor Mj,n and total reflection coefficient Rj,n can be obtained as
For the cylindrical multilayered structure under investigation in this paper, the acoustic field calculation by other methods may lead to singular values when the material parameters of adjacent layers are of substantial difference and the number of the layers becomes large. The proposed method in this paper can solve such problems more efficiently and accurately. Furthermore, the total scattering cross section of the structure is easy to obtain with this method, providing convenience in evaluating the performance of the invisibility cloak and the possibility of a new design approach.
Here, the capability of the multilayered shell structure for reducing the acoustic scattering from a cylindrical sensor system in the 2D problem is demonstrated by using the method above. The background medium is chosen to be water with ρ0 = 998 kg/m3 and κ0 = 2.19 GPa. The material parameters of the SNM are ρα = ρ0, κα = − κ0, and ρβ = − ρ0, κβ = κ0. The whole sensor system is regarded as being a scattering object with radius a′ = 1 cm, in which the effective mass density and the bulk modulus are ρs = ρ0 and κs = (4/9)κ0, respectively. The ratio between the outer radius b and the inner radius a′ of the annular shell is
In order to quantitatively analyze the concealment effect of the proposed structure on the acoustic sensor, and compare the performances of the cloaking shells with different numbers of layers, the variation of total scattering cross-section with the number of layers increasing (up to N = 100) is computed, and the results are shown in Fig.
Now we consider the case where the equivalent density of the scattering object (bare sensor) is not equal to that of the background medium. According to transformation acoustic theory, the phase matching and impedance matching are established at the same time only when the equivalent density of the scatter is equal to that of the surrounding environment. However, as is well known, the phase matching is more important in those structures. Therefore, on the premise of guaranteeing the phase matching, the impedance matching can be sacrificed to some extent to adapt to the unequal density between the scatter and background medium. So we multiply both the density ρs and bulk modulus κs by a robust factor ξ. It can be seen from Fig.
Further, we calculate the acoustic field distribution when ξ values are 0.5 and 2 as shown in Fig.
According to the previous derivation process, it can be found that the selection of single negative material has no bearing on the parameters of the background medium and the scatter. Here, we investigate the scatterings of the multilayer structure under different material parameters of SNM. We consider two cases: ρα = −3ρ0, κα = 2κ0, ρβ = 3ρ0, κβ = − 2κ0, and ρα = − 0.5ρ0, κα = 0.3κ0, ρβ = − 0.5ρ0, κβ = − 0.3κ0. The corresponding acoustic field distributions are shown in Fig.
As we can see from Fig.
Since single negative material achieved by metamaterial structure inevitably produces dissipation, it is necessary to study the scattering of the multilayer structure with dissipation in SNM. So we add a positive imaginary part to the mass density to describe the dissipation in SNM. Then the mass density used in calculation is ρα = (−1 + id)ρ0, ρβ = (1 + id)ρ0, where d is the dissipation factor. We calculate the variation of the total scattering cross section when the dissipation coefficient changes, and the results are shown in Fig.
The corresponding acoustic pressure distributions are shown in Fig.
The previous cloaking devices are designed according to the superlens cloaking theory, and analyzed with the scattering method. In the following we will show the possibility of the invisibility cloak design with the scattering method. As shown above, the scattering cross section is easy to obtain with this scattering method, by searching, in the realizable acoustic parameter space, for the optimum value of the cloaking shell that can make the sensor undetectable to the signals, and the scattering-cancellation-type cloak can be obtained. It is assumed that the sensor can be homogenized as a soft effective medium with acoustic parameters of 0.13ρ0, 0.15κ0, with radius of 0.5 cm. For the case without the cloak, the scattering field of this sensor excited by a plane wave source or cylindrical wave source with a wave length of 1 cm are shown in Figs.
For simplicity, we only use two layers with the same thickness, i.e., 0.2 cm for cloaking, in which case the superlens cloaking method does not work effectively as shown in Fig.
In this paper, the recursive numerical method is employed to study the cylindrical multilayered system and reduces the scattering from an acoustic sensor while allowing it to receive external information. The acoustic pressure field is obviously disturbed by the scattering of the bare sensor system, and the scattering can be cancelled when the sensor system is wrapped in a 10-SNM layer shell. We calculate the acoustic pressure fields in various cases. First, we analyze the scattering characteristics of the cloaking shell with different numbers of bilayers, the total scattering cross section of the structure decreases as the number of bilayers increases. Then we consider the case where the equivalent density of the bare sensor is not equal to the background medium. As expected, the scattered field can still be suppressed by the multilayered shell in this case. Further, we investigate the effect of the SNM parameter on the scattering field. When none of the absolute values of the material parameters of SNM is strictly equal to its counterpart of the background medium, the scattering of the acoustic pressure field emerges, but the scattering can still be suppressed by the multilayer shell. Finally we find that the backscattering nearly approaches to zero when acoustic dissipation induced by SNM is not very strong and some scatterings happen on the other side of the structure. The amplitude of scattering wave increases as the dissipation factor increases. Last but not least, the recursive method is efficient in calculating the scattering cross section and used to design a scattering-cancellation-type invisibility cloak.
1 | |
2 | |
3 | |
4 | |
5 | |
6 | |
7 | |
8 | |
9 | |
10 | |
11 | |
12 | |
13 | |
14 | |
15 | |
16 | |
17 | |
18 | |
19 | |
20 | |
21 | |
22 | |
23 |